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ABSTRACT 

We present the results of 2. 5 -dimensional resistive magnetohydrodynamic (MHD) simulations of the magnetic 
interaction between a weakly magnetized neutron star and its accretion disk. General relativistic effects are 
simulated by using the pseudo-Newtonian potential. We find that well-coUimated jets traveling along the 
rotation axis of the disk are formed by the following mechanism: (1) The magnetic loops connecting the 
neutron star and the disk are twisted due to the differential rotation between the neutron star and the disk. 
(2) Twist injection from the disk initiates expansion of the loop. (3) The expanding magnetic loops create a 
magnetic tower in which accelerated disk material travel as collimated bipolar jets. The propagation speed of 
the working surface of the jet is the order of 10% of the speed of light 0.1c). (4) Magnetic reconnection 
taking place inside the expanding magnetic loops injects hot bubbles intermittently into the magnetic tower. 
The ejection speed of the bubble is the order of the local Alfven speed of the launching point and ^ 0.2c in our 
simulations. (5) The hot bubbles moving inside the tower catch up with the working surface of the jet. High 
energy electrons created by the magnetic reconnection are a plausible source of radio emission. Our model can 
explain the formation process of a narrow jet from a weakly magnetized < 10^ gauss) neutron star and 
the correlation between radio flares of the core and of the lobe observed in Sco X-1. 

Subject headings: accretion, accretion disks — relativity — MHD — stars: neutron — ISM: jets and outflows 



1. INTRODUCTION 

The mechanism of jet formation is one of the most impor- 
tant subjects of research in astrophysics. Blandford & Payne 
(1982) studied magneto-centrifugal acceleration along a mag- 
netic field line threading an accretion disk. Uchida & Shi- 
bata (1985) and Shibata & Uchida (1986) showed by time- 
dependent magnetohydrodynamical (MHD) simulations that 
bipolar jets are produced from an accretion disk threaded by 
open magnetic field lines. 

Observational evidence indicate the presence of semi- 
relativistic jets in some galactic X-ray binaries (XRBs) such 
as SS433 (Hjellming & Johnston 1981), Cygnus X-3 (Marti 
et al. 2000, Mioduszewski et al. 2001), Sco X-1 (Fomalont 
etal. 2001a,b), and Circinus X-1 (Fender ef aZ. 1998). The 
speed of jets in such XRBs indicates that the jets are launched 
close to the last stable orbit of the accretion disk. Therefore, 
we should include general relativistic effects around a com- 
pact object. In addition, the magnetic field of the compact 
object may play an essential role in the dynamics of the jets 
and the accretion flow. Recent observations of XRBs by the 
Very Large Baseline Interferometer (VLBI) revealed their en- 
ergetic activities such as radio flares and radio jets. VLBI 
discovered the motions and variabilities of radio lobes in Sco 
X-1 (Fomalont et al. 2001a,b). According to their observa- 
tions, radio lobes are ejected from the core of Sco X-1. After 
the radio flare of the core, the intensity of the advancing ra- 
dio lobe increases. The correlation between the core flare and 
the lobe flare in Sco X-1 indicates that explosive events at the 
core transports energy to the lobe. 

A plausible model of advancing lobes is the working sur- 
face of jets ejected from the core. However, the existence 
of large scale open magnetic fields assumed in the magneto- 
centrifugal model of jet formation is not evident in XRBs. 
Instead, we should take into account the magnetic interaction 
between the neutron star and its accretion disk. 

Lynden-Bell & Boily (1994) studied the evolution of force- 
free magnetic loops anchored to the star and the disk. They 



obtained self-similar solutions for the evolution of the mag- 
netic loops. They found that the loops are unstable against 
the twist injection from the rotating disk and that the loops 
expand along a direction of 60 degrees from the rotation axis 
of the disk. Lovelace et al. (1995) pointed out that the dipole 
magnetic field of the star deforms itself into an open magnetic 
field due to the differential rotation between the star and the 
disk. Subsequently, Lynden-Bell (1996; hereafter referred to 
as LB96) showed that a cylindrical magnetic tower is formed 
when the disk is surrounded by external plasma with finite 
pressure. 

Hayashi, Shibata, and Matsumoto (1996; hereafter referred 
to as HSM96) carried out the first numerical calculations of 
magnetic interaction between a protostar and its accretion 
disk. By assuming Newtonian gravity and neglecting the ro- 
tation of the star, they showed that the magnetic interaction 
can explain the X-ray flares and outflows observed in pro- 
tostars. Later, Goodson, Bohm, and Winglee; Goodson & 
Winglee (1999); hereafter referred to as GBW99 and GW99, 
respectively, showed that the recurring magnetic reconnection 
creates periodic outflow along the rotation axis of the disk. 
However, these works have not clearly shown the formation 
of a magnetic tower 

In this paper, we numerically studied the formation process 
of magnetic towers around a weakly magnetized neutron star 
In §2 we present basic equations and models. Numerical re- 
sults are shown in §3. §4 is devoted for discussion and con- 
clusion. 

2. SIMULATION MODEL 

We solved resistive magnetohydrodynamic equations in 
cylindrical coordinates. For normalization, we use the 
Schwarzschild radius rj = 2GM/c^ = 1 and the speed of light 
c = 1 . General relativistic effects in the innermost region of 
the accretion disk around the neutron star are taken into ac- 
count by using the pseudo-Newtonian potential (Paczyrisky & 
Wiita 1980) i' = -l /[2{R-1)] where R = y/?^. The basic 
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equations in conservative form are as follows: 

fw.(p.) = o 
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where e = pv^/2 + p/ij- 1) is the total energy of the gas 
(7 = 4/3), E = -V X B + r]J and J = (V X B)/47r. The other 
symbols have their usual meanings. 

The resistivity rj is assumed to be uniform in the en- 
tire computational domain and parametrized by the magnetic 
Reynolds number = crs/r] = 1000, where we take the char- 
acteristic speed and length as c and r^, respectively. The re- 
sistivity is taken to be small enough to ensure the coupling 
between plasma and magnetic fields. The diffusion time scale 
of magnetic fields = A/rj = (A/crs)Rm (A is the area of 
the diffusion region) is much longer than the dynamical time 
scale except in the localized current sheet. In this paper, we 
assumed a uniform resistivity rather than an anomalous resis- 
tivity, which incorporates plasma instabilities due to kinetic 
effects, adopted in Hayashi et al. (1996). As we show later, 
essential results of Hayashi et al.'s simulation such as a topo- 
logical change of magnetic field lines and ejection of plas- 
moids are reproduced by numerical simulations with small 
uniform resistivity. 

The basic equations are solved by the 2-D axisymmetric 
MHD code (Matsumoto et al. 1996; Hayashi et al. 1996) 
based on a modified Lax-Wendroff scheme with artificial vis- 
cosity. 

The initial rotating torus is obtained by assuming a poly- 
tropic equation of state p = Kp^^^l" where n = 3 and the 
angular momentum distribution is I = loir/ro)" where Iq = 

^Q^^/[(fo~ 1)V2] and a = 0.35. The subscript signifies the 
values at the density maximum of the torus ro = 13 and the 
density is normalized at this point. The initial density distri- 
bution of the torus is given by: 
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where po, ip = -ip + {I / r)^ / [2{1 -a)], and Vs,o are the maximum 
density of the torus, the effective potential, and the sound 
speed of the torus, respectively. The thermal energy of the 
torus is parametrized by Eth.o = o/TlV'ol- 
Outside the torus, we assume a non-rotating, spherical, and 
isothermal hot corona in hydrostatic equilibrium. The density 
distribution of the corona is: 



Pc = p* exp 
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The subscript * signifies the values at the surface of the neu- 
tron star = 2.8. The thermal energy of the corona is 
parametrized by isth, * = ^ /7 1 ■0* | ■ 



The initial magnetic field of the neutron star is the dipole 
magnetic field described by the toroidal component of the vec- 
tor potential A ^ = -m^r/R . The strength of the magnetic field 
is defined by the magnetic dipole moment: 

^ / lirp., V?,* 

parametrized by the value of the plasma /3 defined as the ratio 
of gas pressure of the corona to magnetic pressure. 

Initial conditions are determined by choosing the non- 
dimensional parameters p*/po = 0.006, Eth.o = 0.006, E±,* = 
0.21, and /3* = 2. We use 600 x 600 non-uniform meshes. The 
grid size is uniform ( Ar = Az = 0.025) in the region < r < 10 
and < z < 10. Otherwise, it increases with r and z where 
r > 10orz> 10. The size of the entire computational box is 
< r < 100 and < z < 100. At the inner boundary R = R^., 
we imposed an absorbing boundary condition. The deviation 
from the initial values are artificially decreased in the damp- 
ing layer between R = R^, and R = R^ + AR, where AR = 0.2. 
We imposed a symmetric boundary condition at the equato- 
rial plane and a reflecting boundary condition at the rotation 
axis. The outer boundaries are free boundaries where waves 
can transmit. 

3. SIMULATION RESULTS 

Figure[l]shows the time evolution of poloidal magnetic field 
lines (solid lines) and the distribution ofT = 'yP/ p (color con- 
tours). The arrows show the velocity vectors. In the upper 
panels, the time interval between each panel is about one rota- 
tion period r„„ ^ 35Qrs/c at r = ro. The bottom panels enlarge 
the innermost region. 

As the magnetic field lines connecting the neutron star and 
the disk are twisted due to the rotation of the disk, they be- 
gin to inflate in the direction 30-60 degrees from the rotation 
axis. The magnetic fields cease to splay out when the mag- 
netic pressure balances with the ambient gas pressure. After- 
wards, the expanding magnetic loops form a cylindrical tower 
of helical magnetic fields whose height increases with time. In 
previous MHD simulations of disk-star magnetic interaction 
(e.g., HSM96; GBW99 and GW99; Ustyugova et al. 2000), 
the tower structure was not so prominent because the ambient 
gas pressure was too low to confine the magnetic tower inside 
the computational box. 

Numerical results indicate that the top interface between 
the tower and the external medium propagates with a speed 
~ 0.1c. This propagation speed is consistent with that of the 
theory of magnetic towers (LB96). Disk materials are acceler- 
ated and launched along the magnetic tower. Inside the mag- 
netic tower, current sheets are formed because magnetic field 
lines extend upward from the disk and then go downward to 
the star. Magnetic reconnection taking place in the current 
sheet injects hot plasmoids intermittently into the magnetic 
tower. In figure [^, symbols 'B', 'C, and 'D' denote such 
plasmoids created by magnetic reconnection. Region 'A' is 
the hot region between the top of the magnetic tower and the 
jet terminal shock. Beside this working surface, the jet flow 
changes its direction away from the rotation axis and creates 
backflows. The radial interface between the magnetic tower 
and the ambient matter is stable for the growth of the Kelvin- 
Helmholtz instability because magnetic fields stabilize the in- 
stability. 

Figure |2t shows the time evolution of the distribution of 
toroidal magnetic fields. Figure 2b visualizes the three- 
dimensional structure of magnetic field lines at r=739.5. The 
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Fig. 1. — Formation of a magnetic tower, working surface 'A', and intermittent ejection of hot bubbles 'B', 'C, and 'D' and their time evolution inside the 
magnetic tower (upper panels). The color contour shows the temperature distribution and arrows show the velocity vectors. The disk material is launched from 
the inner part of the accretion disk (bottom-left panel) to a direction of 60 degrees from the rotation axis and coUimated into the rotation axis. The initial outflow 
velocity is about 20% of the speed of Ught. Hot bubbles are created inside the magnetic loops anchored to the disk and the star (bottom-right panel.) 



magnetic field lines are strongly twisted due to the rotation 
of the disk. Note the similarity between this figure and the 
analytical model of magnetic towers (figure 2 of LB96.) 

Figure |3^ shows the trajectories of the working surface of 
the magnetic tower 'A' and the bubbles 'B', 'C, and 'D'. The 
plasmoids are ejected intermittently with interval of the ro- 
tation period around 5- lOr,. The maximum speed of the 
plasmoid moving inside the magnetic tower is about 0.2c. 
The plasmoid 'B' catches up with the working surface 'A' 
at t ^ 800 and merges after t ~ 900. 



Figure Isj) shows the distribution of velocity, density, tem- 
perature, gas pressure, magnetic pressure, and the verti- 
cal energy flux [kinetic energy flux (pv^/2)v,, enthalpy flux 
llP/h- ^)]vz^ the Poynting flux (E x B)-/4t:] at r = 10 (in- 
side the magnetic tower) and r = 30 (outside the magnetic 
tower). The radius of the tower is about 20. The total pressure 
of the tower slightly exceeds the ambient gas pressure. It is 
clear that the Poynting flux dominates other fluxes inside the 
jets. 

The density distribution has a minimum in the working sur- 
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Fig. 2. — (a) Magnetic field configuration of tlie magnetic tower in r-i plane (top panels). Solid curves show poloidal magnetic field lines projected onto 
this plane. The color contour shows the strength of toroidal magnetic fields and arrows show velocity vectors. The alphabetical symbols indicate the individual 
magnetic islands, (b) The bottom panel shows the 3-D image of the magnetic field lines (solid lines) and density distribution in the equatorial plane (color 
contour). 



face 'A'. It corresponds to the region between the contact dis- 
continuity (CD) and the jet terminal shock (JTS). We should 
remark that the bow shock ahead of the working surface is 
not visible because the sound speed of the ambient plasma 
(vs „ ~ 0.4) is larger than the propagation speed of the jet. 

4. DISCUSSION 

In this letter, we have demonstrated that highly collimated 
magnetic towers are formed along the rotation axis of the ac- 
cretion disk surrounding a weakly magnetized neutron star. 
We also found that hot plasmoids created by the intermittent 
magnetic reconnections are injected into the magnetic tower 
(Figure 0. The magnetic towers are confined to the radius 
where the magnetic pressure of the expanding magnetic loops 



is comparable to the gas pressure of the ambient plasma. 

Recently, Poynting jet, in which the energy and angular 
momentum are carried predominantly by the electromagnetic 
field, has been studied numerically by Romano va etal. ( 1 998) 
and Ustyugova et al. (2000). The magnetic tower jet is a 
Poynting jet, because the Poynting flux dominates the energy 
flux of the bulk flow inside the jet in our simulation (Figure 
^p). The electromagnetic extraction of angular momentum 
from the disk drives the accretion of the disk material. It 
is interesting to compare the amount of angular momentum 
extracted electromagnetically by the jet and the angular mo- 
mentum transported inside the disk due to the magnetic tur- 
bulence driven by the magneto-rotational instability (MRJ). 
We have to carry out three-dimensional MHD simulations to 
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Fig. 3. — (a) The trajectories of the working surface 'A' and the bubbles 'B', 'C, and 'D' in (t,z) plane. The shaded region indicates the working surface 'A'. 
At f = 900, the plasmoid 'B' colhdes with 'A and merges. The speed of the bubbles is about 10% of the speed of light. The dotted Une indicates 10% of the speed 
of light, (b) The snapshot of velocity, density, temperature, gas pressure (thin), magnetic pressure (thick), and distribution of energy flux inside the magnetic 
tower are shown from top to bottom. Solid and dashed curves indicate the value at r = 10 and r — 30, respectively. A contact discontinuity (CD) locates at z = 67 
and a jet terminal shock (JTS) is placed at z = 52. 



determine the MRI driven angular momentum transport rate 
self-consistently. We intend to report the results of such sim- 
ulations in subsequent papers. Our simulation demonstrated 
the formation mechanism of a semi-relativistic Poynting jet 
around a neutron star 

We can identify the jet terminal shock (JTS) behind the 
working surface of the jet. Our numerical results show that 
a hot, low density region is created between the JTS and the 
discontinuity (CD) separating the jet material and the ambi- 
ent plasma. The structure of the working surface of the jet is 
qualitatively consistent with those in previous MHD simula- 
tions of jet propagation assuming force-free magnetic fields 
initially uniform in the z-direction (e.g., Todo et al. 1992), al- 
though the detailed structure is different because our numeri- 
cal simulations started with a dipole magnetic field. 

We found that hot plasmoids are injected intermittently into 
the magnetic tower. The injection speed of the blobs is the 
order of the local Alfven speed of the reconnection region 
(v^ ^ 0.2c), which is comparable to the rotation speed of the 
innermost region of the disk. Since this speed is faster than 
the propagation speed of the working surface of the jet, the 
blobs catch up with the working surface and release magnetic 
energy. Magnetic reconnection taking place in the core and 
in the working surface could generate high energy electrons 
which emit synchrotron radiation. 

Let us discuss the correlation between the radio flare of the 
core and the lobe observed in Sco X-1. Magnetic reconnec- 
tion taking place in the core will eject hot plasmoids, which 
subsequently release magnetic energy again at the working 
surface of the magnetic tower. This event can be observed 
as the lobe flares correlated with the core flares. In our sim- 
ulations, however, the maximum speed of the hot plasmoids 
(~ 0.2c) is smaller than the speed estimated from the time 



lag between the core flare and the lobe flare. This discrepancy 
may be resolved if the Alfven speed inside the magnetic tower 
is faster than that attained in our simulation. This happens 
when the plasma density inside the magnetic tower is lower 
than the ambient plasma. Plasma heating due to magnetic 
reconnection helps increasing the Alfven speed and ejection 
speed of hot bubbles. 

We should point out the limitation of using non-relativistic 
MHD. Although we take into account general relativistic ef- 
fects by using a pseudo-Newtonian potential, we solved New- 
tonian MHD equations by neglecting special relativistic cor- 
rections. By extending our work to fully general relativis- 
tic MHD, we will be able to simulate models having Alfven 
speed closer to the light speed. 

The work presented here is the first numerical simulation 
which clearly shows the formation process of the axisymmet- 
ric magnetic tower and semi-relativistic jets. We expect that 
semi-relativistic jets observed in some XRBs are driven by 
the magnetic interaction between the dipole magnetic field of 
a weakly magnetized neutron star and its accretion disk. 3-D 
simulations are necessary to study the stability of magnetic 
towers against non-axisymmetric perturbations. We would 
like to report the results of such simulations in the near fu- 
ture. 
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